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ABSTRACT 

We reconstruct the primordial spectrum of the curvature perturbation, P(k), 
from the observational data of the Wilkinson Microwave Anisotropy Probe 
(WMAP) by the cosmic inversion method developed recently. In contrast to 
conventional parameter-fitting methods, our method can potentially reproduce 
small features in P(k) with good accuracy. As a result, we obtain a compli- 
cated oscillatory P(k). We confirm that this reconstructed P(k) recovers the 
WMAP angular power spectrum with resolution up to A£ ~ 5. Similar oscilla- 
tory features are found, however, in simulations using artificial cosmic microwave 
background data generated from a scale-invariant P(k) with random errors that 
mimic observation. In order to examine the statistical significance of the nontriv- 
ial features, including the oscillatory behaviors, therefore, we consider a method 
to quantify the deviation from scale-invariance and apply it to the P(k) recon- 
structed from the WMAP data. We find that there are possible deviations from 
scale-invariance around k ~ 1.5 x 1CT 2 and 2.6 x 10~ 2 Mpc~ . 

Subject headings: cosmic microwave background — cosmology: theory 



1. Introduction 

The Wilkinson Microwave Anisotropy Probe (WMAP) satellite has brought us inter- 
esting information about our universe (Bennett et al. 2003). From the remarkably precise 
observation of the temperature fluctuations and the polarization of the cosmic microwave 
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background (CMB), one can obtain not only accurate values of the global cosmological pa- 
rameters, but also invaluable information of the properties on the primordial fluctuations 
(Bennett et al. 2003; Spergel et al. 2003; Peiris et al. 2003; Komatsu et al. 2003). Although 
their results as a whole support the standard ACDM (cold dark matter) model with Gaus- 
sian, adiabatic, and scale-invariant primordial fluctuations, some features that cannot be 
explained by the standard model have also been pointed out, such as (1) lack of power on 
large scales, (2) running of the spectral index n s from n s > 1 on larger scales to n s < 1 on 
smaller scales, and (3) oscillatory behaviors of the power spectrum on intermediate scales. 

In fact, feature 1 was known already with COBE data (Bennett et al. 1996), and a 
possible explanation was proposed (Yokoyama 1999). There have been many new proposals 
(Bi et al. 2003; de Deo et al. 2003; Efstathiou 2003a; Uzan et al. 2004; Contaldi et al. 2003; 
Cline et al. 2003; Feng and Zhang 2003; Kawasaki and Takahashi 2003; Dvali and Kachru 
2003; Piao et al. 2003), as well as many arguments about its statistical significance (Spergel et 
al. 2003; Tegmark et al. 2003; de Oliveira-Costa et al. 2003; Gaztahaga et al. 2003; Niarchou 
et al. 2004; Efstathiou 2003b, 2004). A number of inflation models that can account for 
feature 2 have also been proposed recently (Feng et al. 2003; Kyae and Shafi 2003; Kawasaki 
et al. 2003; Huang and Li 2003; Chung et al. 2003; Yamaguchi and Yokoyama 2003; Dvali 
and Kachru 2003). 

On the other hand, feature 3, namely, possible oscillatory behaviors around a simple 
power-law spectrum, is more difficult to quantify (Peiris et al. 2003). Although several 
attempts to reconstruct the primordial spectrum were made by combining WMAP data 
with other independent observational data, such as the Two-Degree Field Galaxy Redshift 
Survey (2dFGRS) and Lya forest data (Wang et al. 1999a,b; Hannestad 2001, 2003; Tegmark 
and Zaldarriaga 2002; Wang and Mathews 2002; Bridle et al. 2003; Mukherjee and Wang 
2003a,b,c), they all employed the binning, wavelet band powers, or direct wavelet expansion 
method to the data. These methods cannot detect possible oscillatory behaviors if their 
scale is smaller than the binning scale. It is preferable to use a method that can restore the 
primordial spectrum as a continuous function without any ad hoc filtering scale to investigate 
detailed features such as 3. Fortunately, such a new method has been proposed in Matsumiya 
et al. (2002) and Matsumiya et al. (2003), with the name cosmic inversion method, and their 
test calculations using mock data without observational errors have shown that this method 
can reproduce possible small dips and peaks off a scale-invariant spectrum quite well, in 
the spatially flat universe with the adiabatic initial condition. In this paper, we attempt 
to reconstruct P(k) from the WMAP data using this method. To apply it to the actual 
data, we consider the observational errors by Monte Carlo simulations and use a modified 
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CMBFAST 1 code that adopts much finer resolution than the original one in both k and £, so 
that we can compute the fine structure of angular power spectra accurately. We mention that 
recently another new method for the reconstruction of P(k) has been proposed by Shafieloo 
and Souradeep (2003). 

This paper is organized as follows. In Sec. 2, we review the cosmic inversion method 
and describe our analysis method using the WMAP data. In Sec. 3, we apply our method to 
the WMAP data and discuss the results. In particular, we investigate statistical significance 
of deviations from scale-invariance. Finally, we present our conclusion in Sec. 4. 



2. Inversion Method 

First, we briefly review the cosmic inversion method proposed by Matsumiya et al. 
(2002, 2003). We assume a spatially flat universe with an adiabatic initial condition, both of 
which are generic predictions of standard inflation and have been supported by the WMAP 
data. In the end, we obtain a first-order differential equation for the primordial spectrum 
P(k). 

The CMB anisotropy is quantified by the angular correlation function defined as 

C(0) = (e( ni )e(n 2 )) = -^CfPtXcosd), costf = ni • n 2 , (1) 

where 0(n) is the temperature fluctuation in the direction n. We decompose the Fourier 
components of the temperature fluctuations 0(7/, k) into multipole moments, 

oo 

Q( V ,k,n) = J2HY®t(v,k)Pi(ri, (2) 

where /i = k- n, k is the comoving wavenumber, and r/ is the conformal time, with its present 
value being r] . Using 0^(r?, k), the angular power spectrum is expressed as 

2£+l^ 1 f°°dk k 3 {\&e(vo,k)\ 2 } 
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The Boltzmann equation for 0(r/, k) can be transformed into the following integral form (Hu 
and Sugiyama 1995): 

(e + tf)(77o,fc,AO = J rf»7 1 [0 O + ^ -ifJ.V h ]V(ri) + (V - $) e -^>) je^-^ ( 4 ) 



1 See http://www.cmbfast.org/ 
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where the overdot denotes the derivative with respect to the conformal time. Here, ^ and 
$ are the Newton potential and the spatial curvature perturbation in the Newton gauge, 
respectively (Kodama and Sasaki 1984), and 



-t(tj) 



V(r)) = te- T{r <\ r(r])= / idr], (5) 

are the visibility function and the optical depth for Thomson scattering, respectively. In the 
limit that the thickness of the last scattering surface (LSS) is negligible, we have V{rf) m 
5(i] — 77*), e~ T ^ ~ 6{r] — 77*), where 77* is the recombination time when the visibility function 
is maximum (Hu and Sugiyama 1995). Taking the thickness of the LSS into account, we 
have a better approximation for Eq. (4) as 

(e + *)(7to,A;,f/)« / rfr/|[e o + ^-^0i]v(r/) + (^-<i>)e- r(r ' ) U-^ = e app + *, (6) 

J ??*start 

where d = i]q — i]* is the conformal distance from the present to the LSS and restart and r?* en d 
are the times when the recombination starts and ends, respectively. Here we introduce the 
transfer functions, f{k) and g(k), defined by 

d V [(Q + ^)( V ,k)V( V ) + ^ -^>)( V ,k)e- T ^\ = /(A;)$(0,k), (7) 

' ?7*start 

/ d7 7 e 1 (7 7 ,A;)V(77) ee <?(£;)<t>(0,k). (8) 

??*start 

We can calculate /(£;) and g{k) numerically; they depend only on the cosmological param- 
eters, for we are assuming that only adiabatic fluctuations are present. Then, we find the 
approximated multipole moments as 



ef p (77o, k) = (21 + 1) [f(k)j e (kd) + g(k)j' e (kd)} $(0, k), (9) 
and the approximated angular correlation function as 

C-(r) = £ ^Q app P, (l - , (10) 

where C app is obtained by putting Eq. (9) into Eq. (3), r is defined as r = 2dsin(6/2) on 
the LSS, and £ min and £ max are lower and upper bounds on £ due to the limitation of the 
observation. In the small-scale limit r <C d, using the Fourier sine formula, we obtain a first- 
order differential equation for the primordial power spectrum of the curvature perturbation, 

P(fc) = <|$(0,k)| 2 >, 

-k 2 f(k)P'(k) + [-2k 2 f{k)f\k) + kg\k)] P{k) 
r°° 1 f) 

= 4ir dr-— {r 3 C app (r)} sin kr = S(k). (11) 
Jo rdr 
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Since f(k) and g(k) are oscillatory functions around zero, we can find values of P{k) at the 
zero points of f(k) as 

P{ks) = i^th for f{ks)=0 > (12) 

assuming that P'(k) is finite at the singularities, k = k s . If the cosmological parameters and 
the angular power spectrum are given, we can solve Eq. (11) as a boundary value problem 
between the singularities. 

However, because Eq. (11) is derived by adopting the approximation in Eq. (6), C^ pp 
is different from the exact angular spectrum C| x for the same initial spectrum. The errors 
caused by the approximation, or the relative differences between C^ pp and C| x , are as large as 
about 30%. Thus, we should not use the observed power spectrum C obs directly in Eq. (10). 
Instead, we must find the Cf v that would be obtained for the actual P(k). Of course, this 
is impossible in the rigorous sense, since it is the actual P{k) that we are to reconstruct. 
It turns out to be possible, however, with accuracy high enough for our present purpose, 
because we find that the ratio, 

fen 

be = T^p, (13) 

is almost independent of P(k) (Matsumiya et al. 2003). Using this fact, we first calculate 
the ratio, bf^ = Cj x ^ / C^ pp< "°\ for a known fiducial initial spectrum P(°\k) such as the 
scale-invariant one. Then, inserting Cf 38 /bf\ which is much closer to the actual Cf p , into 
the source term of Eq. (11), we may solve for P(k) with good accuracy. 

In practice, we cannot take the upper bound of the integration in the right-hand side of 
Eq. (11) to be infinite. The integrand in Eq. (11) is oscillating with its amplitude increasing 
with r. We therefore introduce a cutoff scale r cut . However, this inevitably introduces a 
smoothing scale to our method. As the cutoff scale is made larger, the rapid oscillations of 
the integrand with increasing amplitude become numerically uncontrollable. On the other 
hand, if the cutoff scale is made smaller, the resolution in the /c-space becomes worse as 
A A; ~ n/r cnt . In the actual calculations, to maintain as good a numerical accuracy as 
possible and to obtain simultaneously as fine a resolution in /c-space as possible, we convolve 
an exponentially decreasing function into the integration of the Fourier sine transform with 
the optimized cutoff scale of r cut ~ 0.5d, corresponding to 9 ~ 30°, or Al ~ 6. Thus, the 
resolution of the Fourier sine transform is limited to Akd ~ 6. We could not adopt a finer 
resolution because of the numerical instability. Nevertheless, the resolution of our method is 
much finer than that of any other method. Note that there is an absolute theoretical limit 
^cut ^ 2d, or Akd > 1.5, because of the finite size of the LSS sphere. 
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In this paper, we also account for the effect of observational errors on the reconstructed 
P(k). First, for each £ we draw a random number from a Gaussian distribution whose mean 
value is equal to a central value of the observed and variance is given by the diagonal 
term of the covariance matrix, and then we reconstruct P(k) from each simulated data 
set. In fact, each Ci is weakly correlated with other multipoles and follows a x 2 distribution. 
However, our procedure is valid because the correlation with other multipoles is weak enough 
and the x 2 distribution is practically identical to the Gaussian distribution for sufficiently 
large £ where we analyze. We estimate the mean value and the variance of the reconstructed 
P(k) at each k for 1000 realizations. We find that the intrinsic errors caused by our inversion 
method itself, whose magnitude is estimated by inverting Ci spectra calculated from artificial 
P{k) spectra without observational errors, are much smaller than the observational errors of 
WMAP, except around the singularities where the numerical errors are amplified. 

To calculate C£ x , we use a modified CMBFAST code with much finer resolutions than 
the original one in both k and I. We limit Ce in the range 20 < £ < 700 in order not 
to use data that have large observational errors due to the cosmic variance at small £ and 
the detector noise at large £. We adopt the fiducial initial spectrum, P^°\k), as the scale- 
invariant spectrum, and the fiducial cosmological parameter set as h = 0.72, Q b = 0.047, 
Q\ = 0.71, Q m = 0.29, and r = 0.17, which are the best-fit values to the WMAP data for the 
scale-invariant spectrum, k 3 P(k) = A. In this case, the positions of the singularities given by 
Eq. (12) are kd ~ 70, 430, 680, where d ~ 1.34 x 10 4 Mpc. Because the reconstructed P(k) 
around the singularities has large numerical errors that are amplified by the observational 
errors, we can obtain P(k) with good accuracy in the limited range 120 < kd < 380 or 9.0 x 
lO^Mpc^ 1 < k < 2.8xlO~ 2 Mpc~ 1 . Then, we vary some cosmological parameters to examine 
how the shape of the reconstructed P(k) is affected. That is, we also calculate for some 
cases of h = 0.65,0.70,0.75,0.80, tt b = 0.03,0.04,0.05,0.06, and fi A = 0.65,0.70,0.75,0.80, 
respectively. Note that since r affects the shape of the power spectrum only on large scales, 
except for the normalization, we use a fixed value as r = 0.17. 



3. Results and Discussion 

3.1. Reconstruction from Original Data 

We show P(k) reconstructed from the WMAP data for the fiducial cosmological param- 
eter set, namely, h = 0.72, tt b = 0.047, fi A = 0.71, Vt m = 0.29, and r = 0.17 in Fig. 1. We 
can see oscillations whose amplitude is about 20% — 30% of the mean value with frequency 
(Akdy 1 ~ 1/15 — 1/10. To check whether our method works correctly, we recalculate C e 
from the obtained P(k) in the range 120 < kd < 380, assuming scale-invariance outside of 
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this range, and compare it with the observational data of WMAP. As mentioned in Sec. 2, our 
method reconstructs P(k) with a finite resolution, which is caused by the cutoff scale of the 
Fourier sine integral. In Fig. 2, we compare the cases of r cut ~ 0.5c? and ~ 0.3d, which lead 
to A£ ~ Akd ~ 6 or AA; ~ 4.5 x 10 -4 Mpc -1 , and A£ ~ Akd ~ 10 or A A; ~ 7.5 x 10 -4 Mpc -1 , 
respectively. We find that the recalculated Ce spectra agree with the binned WMAP data 
corresponding to the respective smoothing scales as shown in Fig. 3. These agreements are 
quite impressive. We also note that the characteristic frequency of oscillation changes ac- 
cordingly as we vary the cutoff scale, so the observed oscillatory behavior with frequency 
(Akdy 1 ~ 1/15 — 1/10 in Fig. 1 does not necessarily have fundamental meaning. We confirm 
that the resolution becomes better as the cutoff scale is made larger, in agreement with the 
relation Ak ~ 7r/r cut . Of course, we adopt the former resolution, the finest possible scheme 
without numerical instability in subsequent sections. 



3.2. Reconstruction from Binned Data 

Figure 4 demonstrates the dependence of the reconstructed P(k) on the bin size of the 
data. These P(k) spectra are reconstructed from the binned WMAP data whose bin sizes 
are A£ = 10, 20, and 50, by interpolating through the binned data points, respectively. This 
is to see both global and local features in P(k). It is found that the oscillatory feature 
becomes more prominent as the bin size is made smaller, even if the spectrum globally looks 
scale-invariant. We emphasize that such nontrivial features cannot be quantified in P(k) as 
long as conventional parameter-fitting methods are used. 



3.3. Statistical Analysis 

To examine whether the oscillatory features are real, we perform the following sim- 
ulations and compare the results. First, we assume a scale-invariant P(k) and calculate 
theoretical Ce spectra for the fiducial cosmological parameter set. Then, we make artificial 
data of Ce spectra with the same errors as those of the WMAP data by drawing random 
numbers. We use these artificial data to reconstruct P(k), assuming that the cosmological 
parameters are known. The resultant P(k) spectra for some different realizations are shown 
in Fig. 5. We can see that these reconstructed spectra also have oscillatory features whose 
amplitude and frequency are almost the same as the P(k) from the WMAP data. This is 
caused purely by the scatter of the data. In other words, even if P(k) is really scale- invariant, 
it is likely that the reconstructed P(k) looks oscillatory because of the observational errors, 
and thus it is difficult to conclude whether or not there are some significant features in the 
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reconstructed P(k) from the WMAP data. It is necessary to quantify their significance. 

For this purpose, we define the deviation from the scale-invariance in P(k) in the range 
between k\ and k 2 as 



and evaluate their statistical significance as follows. First, we calculate D(ki,k 2 ) for each 
reconstructed P(k) from mock Ci data for the scale-invariant P(k) in the same simulation 
as mentioned above, with the fiducial cosmological parameters. Then, we estimate the 
probability that D(k 1 ,k 2 ) for the simulation exceeds its observed value in the same range 
of k, by counting such events for 10,000 realizations. We show the results of this analysis 
in Table 1. It is found that only 1.41% of the simulations exceed the observed value in the 
range 100 < kd < 400, and especially, only 0.39% in the range 175 < kd < 225 and 0.59% 
in the range 325 < kd < 375, from which we may conclude that there are some possible 
deviations around kd ~ 200 and ~ 350, as we can also see some prominent features in Fig. 1. 
This is consistent with Spergel et al. (2003), who show the contribution to x 2 P er multipole 
in Ce after fitting the ACDM model with the power law P(k) to the data. 



As is seen in Fig. 1, the reconstructed P(k) exhibits severe oscillations, reaching negative 
values locally around the singularity kd ~ 430. One may wonder whether this is due to an 
inappropriate choice of the fiducial cosmological parameters, because if the different cosmo- 
logical parameters from actual values are adopted in the inversion, the reconstructed P(k) 
is distorted from the real one, especially around the singularities as shown by Matsumiya 
et al. (2003). We find, however, that 39.36% of the simulations based on the scale-invariant 
spectrum with the cosmological parameters fixed to the fiducial values result in larger values 
of D(ki,k 2 ) in the range 380 < kd < 430. Hence, it is likely that the severe oscillations 
observed in Fig. 1 are not caused by an inappropriate choice of the cosmological parameters 
but are simply due to large observational errors. Still, it is important to see the dependence 
of the reconstructed P(k) on the cosmological parameters, for the oscillatory features may 
change and we may be able to find a better choice of their values. In the ideal situations with 
sufficiently small observational errors, we may even constrain the cosmological parameters 
by requiring that the reconstructed P(k) be positive definite. 

Figure 6 depicts the dependence of the reconstructed P(k) on the cosmological parame- 
ters. It is found that the oscillatory features remain, but the global amplitude and especially 




(14) 



3.4. Cosmological Parameters 
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the sharpness of the oscillations around the second singularity change. We also find some 
degree of degeneracy between h, and Q\. In Fig. 6, we can see that the increase in h has 
a similar effect as the decrease in f2& or that is, the amplitude of the oscillations increases 
similarly. This degeneracy, however, may be resolved if we could observe fine structures. For 
example, the reconstructed P(k) is practically insensitive to h in the range kd < 250; hence, 
variation in h may be distinguished from that of Qj, or Q/y. 

With the current magnitude of observational errors in Cg, we encounter negative values 
of the reconstructed P(k) for the scale-invariant spectrum. Hence, we need Cg with smaller 
errors to apply this positivity criterion directly. In any case, we should examine P(k) sys- 
tematically for as many parameter sets as possible beyond those shown in Fig. 6, and we 
need a statistical method to quantify the positivity criterion in the presence of observational 
errors. This issue is currently under study. 

At present, we can reconstruct P(k) only in the range 100 < kd < 400. If Ce data 
improve at £ > 700, we will be able to obtain P(k) up to the third singularity, kd ~ 680. 
This additional information on P(k) around the third singularity will be important for 
constraining the cosmological parameters, partly because of reduction in observational errors 
as discussed above and partly because the shape of P(k) around the singularities is sensitive 
to the cosmological parameters, as mentioned in Sec. 3.3. 

4. Conclusion 

We have reconstructed the shape of the primordial spectrum, P(k), from the WMAP 
angular power spectrum data, Cg, by the inversion method proposed by Matsumiya et al. 
(2002, 2003), under the assumptions that the universe is spatially flat and the primordial 
fluctuations are purely adiabatic. First, we have set the cosmological parameters as h — 0.72, 
Q b = 0.047, Q A = 0.71, Q m = 0.29, and r = 0.17. We have obtained an oscillatory P(k). 
The amplitude of oscillations is about 20% — 30% of the mean value, and the frequency is 
about (Akd)^ 1 ~ 1/15 — 1/10, which reflects the resolution of our method. We have also 
confirmed that the reconstructed P(k) has good accuracy in the range 120 < kd < 380, 
or 9.0 x lO^Mpc -1 < k < 2.8 x lO^Mpc" 1 , with a resolution of Akd ~ 5, or A A; ~ 
3.7 x 10~ 4 Mpc _1 . Thus, our inversion method can reconstruct P(k) with a much finer 
resolution than other methods proposed so far, such as the binning, wavelet band powers, 
or direct wavelet expansion method (Wang et al. 1999a,b; Hannestad 2001, 2003; Wang and 
Mathews 2002; Bridle et al. 2003; Mukherjee and Wang 2003a,b,c). 

To examine the statistical significance of possible nontrivial features in the reconstructed 
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spectrum, we have generated 10,000 sets of mock CMB data from the scale-invariant spec- 
trum and reconstructed P(k) from these artificial data. We have found that only 1.41% of 
them have the values of D(k±, A^) larger than the observed value in the range 100 < kd < 400, 
or 7.5 x lO^Mpc" 1 < k < 3.0 x 10- 2 Mpc _1 , only 0.39% in the range 175 < kd < 225, and 
0.59% in the range 325 < kd < 375. From these results we conclude that there are some 
possible deviations from scale-invariance around kd ~ 200, or k ~ 1.5 x 10~ 2 Mpc _1 , and 
kd ~ 350, or k ~ 2.6 x 10~ 2 Mpc _1 . On the other hand, we conclude that the severe os- 
cillation of the reconstructed P{k) around the singularity kd ~ 430, which drives P{k) to 
negative values in some regions, is due not to an inappropriate choice of the cosmological 
parameters but to the large observational errors. This is because we find a high probability 
around the singularity kd ~ 430, or 39.36% in the range 380 < kd < 430. 

There are some issues that we plan to investigate in detail. First, it will be interesting to 
perform the reconstruction for a much wider range of the cosmological parameters and obtain 
systematic constraints on them, as mentioned in Sec. 3.4. Second, it may be necessary to 
apply the D(k\, /C2) test to other smooth forms of P(k), for example, to those reconstructed 
from the binned data instead of the scale-invariant P{k), as discussed in Sec. 3.2, for various 
values of the cosmological parameters. This is because the global shape of the reconstructed 
P{k) changes depending on the cosmological parameters. Third, we should extend our for- 
malism to include the CMB polarization so that it may be applied not only to the WMAP 
data but also to future CMB observations, such as Planck. 2 Inclusion of the CMB polariza- 
tion will constrain P{k) more severely, and if the B-mode polarization is detected, we can 
investigate the tensor mode of the primordial perturbations (Starobinsky 1979; Rubakov et 
al. 1982; Abbott and Wise 1984; Polnarev 1985; Crittenden et al. 1993a,b). 

We would like to thank Eiichiro Komatsu for fruitful discussions. This work was sup- 
ported in part by JSPS Grants- in- Aid for Scientific Research 12640269 (M.S.) and 13640285 
(J.Y.) and by Monbu-Kagakusho Grant-in-Aid for Scientific Research (S) 14102004 (M.S.). 
N.K. is supported by Research Fellowships of JSPS for Young Scientists (04249). 
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Fig. 1. — Primordial spectrum P(k) reconstructed from the WMAP data for h = 0.72, 
fl& = 0.047, Q\ = 0.71, Q m = 0.29, and r = 0.17. The solid curve and the dashed curves 
represent the mean and the la errors, respectively, of the reconstructed P(k), while the dash- 
dotted curves represent the la from the scale-invariance. The horizontal axis kd corresponds 
roughly to £. The singularities lie at kd ~ 70 and 430. Some prominent features are seen 
around kd ~ 200 and 350. 
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Fig. 2. — Dependence on the cutoff scale r cut . Top: Reconstructed P(k) spectra from the 
WMAP data for r cut — 0.5c? and 0.3d. We see that the resolution becomes better as the 
cutoff scale is made larger. Bottom: Recovered Ce spectra from the obtained P(k) spectra. 
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Fig. 3. — Accuracy check of our reconstruction method. We show the cases of r cut ~ 0.5c/, 
which leads to A£ ~ 6 (top), and r cut ~ 0.3d, which leads to A£ ~ 10 (bottom), from the 
relation A£ ~ A/cd ~ ird/r cut . Left: Comparison of the binned WMAP data, C\™ (plus 
signs with error bars), with the angular power spectrum, C\ e (solid curve), recovered from 
the reconstructed P(k) shown in Fig. 2. Right Relative errors, (Cg — C^ c )/C| c . The 
relative errors are small for most of the bins except for those corresponding to the scales of 
the singularities. 
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Fig. 4. — Reconstruction from the binned WMAP data. Top: Binned WMAP data whose 
bin sizes are A£ = 10, 20, and 50; the solid curve represents the fiducial model for the scale- 
invariant spectrum. Bottom: The P{k) spectra reconstructed from the binned WMAP data. 
As the bin size is made smaller, we can see more oscillations. 
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Fig. 5. — Primordial spectra P(k) reconstructed from artificial CMB data for four different 
realizations. The original spectrum is taken to be scale invariant, and each realization is 
generated by drawing a random number to each Ci with the same error as the WMAP data, 
assuming that the cosmological parameters are known. The same oscillatory features as 
shown in Fig. 1 are seen. 
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Table 1: Statistical Significance of the Deviations from the Scale Invariance 



range [kid, k^d\ 


[100,150 


[150,200 


[200,250] 


[250,300; 


[300,350 


[350,400] 


probability 


18.20% 


9.72% 


3.08% 


83.31% 


16.52% 


2.23% 
















range [k\d, k^d\ 


[100,400 


[175,225] 


[325,375] 


[380,430] 




probability 


1.41% 


0.39% 


0.59% 


39.36% 





Note. — We show the probabilities that the deviation D(k\, k 2 ) for the simulation exceeds its observed 
value in various ranges [kid, k 2 d\. The simulations are performed for the scale-invariant spectrum with the 
fiducial cosmological parameters, which are assumed to be known. 
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Fig. 6. — Some cases of different cosmological parameter sets. We show the reconstructed 
P(k) from the WMAP data for h = 0.65,0.70,0.75,0.80, tt b = 0.04, fi A = 0.70 (top left 
panel), h = 0.70, Q b = 0.03,0.04,0.05,0.06, Q A = 0.70 (top right panel), and h = 0.70, 
Qb = 0.04, = 0.65,0.70,0.75,0.80 (bottom panel), respectively. The oscillatory features 
are roughly similar to each other throughout this range, but they are amplified around the 
second singularity in some cases. 



